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Abstract 

Construction of astrophysically realistic initial data remains a central problem when modelling 
the merger and eventual coalescence of binary black holes in numerical relativity. The objective of 
this paper is to provide astrophysically realistic freely specifiable initial data for binary black hole 
systems in numerical relativity, which are in agreement with post-Newtonian results. Following 
the approach taken by Blanchet [jj, we propose a particular solution to the time- asymmetric 
constraint equations, which represent a system of two moving black holes, in the form of the 
standard conformal decomposition of the spatial metric and the extrinsic curvature. The solution 
for the spatial metric is given in symmetric tracefree form, as well as in Dirac coordinates. We show 
that the solution differs from the usual post-Newtonian metric up to the 2PN order by a coordinate 
transformation. In addition, the solutions, defined at every point of space, differ at second post- 
Newtonian order from the exact, conformally flat, Bowen-York solution of the constraints. 
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I. INTRODUCTION 



General Context 



Construction of 'astrophysically motivated' initial data is of utmost importance in efforts 
directed towards simulating binary black hole (BBH) mergers. Such simulations are crucial 
for the detection and characterization of gravitational waves by the currently operational 
ground-based (e.g. LIGO, VIRGO, TAMA, GEO600), future advanced ground-based and 
space-based (e.g. LISA) laser interferometers. The sensitivity of the detectors is considerably 
increased by matched filtering of the observed signal with sets of the expected full waveform. 
At present, the entire inspiral and eventual coalescence of the two bodies is viewed in terms 
of three consecutive phases. The first and longest duration phase, the gradual 'adiabatic' 
inspiral, is modelled accurately by post-Newtonian (PN) methods and solved analytically 
;o an accuracy of 3.5PN 1 (in the equations of motion and gravitational flux respectively) 
When the PN approximation is assumed to be no longer valid, numerical 
relativity simulations are used at present to model the second phase of the plunge and 



coalescence of the two bodies 
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19, 



20j. After the 



wo compact bodies have merged, the resulting single object enters the final ringdown phase 



2l| . modelled accurately using linear perturbation techniques of single black holes. 



Recently, significant advances have been made in the numerical simulation of the strong- 
field merger phase 14, 3, (ill Q, [3, Q , in particular, in its dynamical evolution. The 
latter enables the accurate description of the BBH's dynamical evolution and the extraction 
of an emitted gravitational wave signal. In this respect, 'astrophysically realistic' initial data 
with an appropriate gravitational wave content would help to ensure that the final extracted 
signal is as close as desirable to the emitted waveform. Note that the term 'astrophysically 
motivated, or realistic, data' refers to initial data sets that are constructed from spacetimes 
representing astrophysical BBHs, the merger of which current and future generation inter- 
ferometers will be able to detect after several centuries or millennia of their inspiral. As 
noted in Section II Bl initial data sets used in early works have not taken into account the 
information about the inspiral and radiation content of BBH systems prior to merger. 



1 The notation 1PN corresponds to the formal ~ 1/c 2 level in a post-Newtonian expansion with respect to 
the Newtonian acceleration and gravitational flux (where c is the speed of light). 
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Initial data is given traditionally in the standard 3+1 formalism |22j and its construction 
is non-trivial. The problem lies in choosing physically meaningful initial data, amongst 
an infinitely large number of non-physical data choices, representing the BBHs. The most 
widely used methods to construct initial data, such as the Extrinsic Curvature (EC) 2^ 
and Conformal Thin-Sandwich (CTS) 24|, adopt conformal decompositions of both the 
spatial 3-metric and the extrinsic curvature. This results in a set of coupled elliptic partial 
differential equations, which are to be solved numerically under appropriate inner and outer 
boundary conditions, together with a set of freely specifiable parameters, referred to later 
in this paper as 'free data'. 



B. Related Work 

For numerical convenience, early works [25. 26] devoted 1 o construct ing initial dat a. sets for 
BBHs introduce a number of significant simplifying assumptions. These include conformal 
flatness, maximal slicing, together with the Bowen-York extrinsic curvatur e 2 |27| . However, 
with considerable recent advances in spectral methods, as exemplified by ^igi bd. . and 
finite element and finite difference techniques, limitations of numerical approaches are no 
onger seen as a serious obstacle. Moreover, several highly-developed 3D numerical software 

dQy y y y a for ^ black hol e ^ ^ ^ ^ — 

wave extraction. This makes it timely now to go beyond the initial assumptions of conformal 
flatness, maximal slicing and the use of the Bowen-York extrinsic curvature, and to provide 
astrophysically motivated initial data. 

Am ong the more recent works that aim to overcome the limitations of the above assump- 



tions, 1 31 



32 



use the alternative double superposed Kerr-Schi 



u tippioaci 

3, HQ 



d approach, [34] uses the 
use an asymp- 



PN method in ADM Transverse Traceless ( ADMTT) coordinates 

totic matching scheme, and j3, Q,S,S] develop helical-Killing vector /quasi-equilibrium 
initial data. These four approaches are briefly discussed below; the reader is referred to 
and for a more detailed discussion. 

Taking the double superposed Kerr-Schild approach first, proposed initially in jjyj and 



2 This is the analytical solution to the decoupled momentum constraint for a conformally flat spacetimc 
and maximal hypcrsurface, was first given in _27], where it was checked by two surface integrals over a 
sphere at infinity. 
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developed further in [32|, |33J, it is important to note that Kerr spacetimes do not admit 
conformally flat slices. In this approach, the free data are constructed from a linear superpo- 
sition of two single (spinning or boosted) black holes in Kerr-Schild coordinates representing 
the spatial conformal metric. Furthermore, the Kerr-Shild data advantageously reduce into 
stationary Kerr in the specific limit of two widely-separated black holes. The Kerr-Schild 



free data have also been numerically implemented 



32« to obtain full initial data sets 3 and 



are amongst the more extensively studied of all astrophysically motivated free data; see, for 



example, 



42|. However, as this work notes (see Appendix |XJ) , the superposed Kerr-Schild 



initial data disagrees with post-Newtonian results at the second-post-Newtonian order. 



In relation to PN free data, similar to Kerr spacetimes, 



43, 



ainary black hole systems are 



44j . It is widely acknowledged 



also not conformally flat at second post-Newtonian order 
now that post-Newtonian results modelling the gradual inspiral should suggest a more astro- 
physically relevant form of the free data for the last few stable orbits prior to the Innermost 
Circular Orbit (ICO) 4 . This is an important motivation for our work, and has also been 
considered in the PN method in ADMTT coordinates j^j]. 

The PN method in ADMTT coordinates due to Tichy et al. proposes the direct 
adoption of the 2PN metric in ADM coordinates as the conformal spatial metric. The 



scheme outlined in |34| is for a numerically computed conformal factor to 'correct' the initial 
post-Newtonian solution and, thus, to account for higher-order post-Newtonian terms in the 
physical metric, which were initially neglected in the proposed conformal counterpart. The 
PN initial data in 2^ are based on the PN results in the ADMTT gauge. Motivations given 
there for the choice of ADMTT gauge apply equally well to PN free data in the harmonic 

l&uge, which is used in our work. Among these commonalities between the two methods, i.e. 

34j and our work, firstly, both methods easily find expressions for the conformal 3-metric and 
extrinsic curvature without any logarithmic terms. Secondly, to the 1.5PN order ~ O(js), 

;he freely specifiable data sets due to both methods also agree with the puncture approach 

2f| and the choice in a specific hypersurface where the trace of the extrinsic curvature, K, 



3 These represent the solution, obtained numerically, to the coupled elliptic partial differential constraint 

equations using the freely specifiable parameters, or free data. 

4 The ICO is defined by the minimum of the binary's energy function for circular orbits. The definition 
does not include, significantly, radiation reaction terms and, therefore, is not physically meaningful in 
the context of exact radiative solutions. It is well-defined, however, for both numerical 00 and 3PN 
calculations, and allows for comparisons between the two facilitating greater understanding. 
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(otherwise known as the 'mean curvature') vanishes. Finally, the two free data sets do not 
include, for the time being at least, any radiation content. 

In the asymptotic matching scheme, following the initial work by Alvi 35], Yunes et al. 



36 



propose an approximate metric for corotating BBHs spacetimes. Advantageously, 
this includes both 1PN effects, and, unlike in post-Newtonian free data sets due to Tichy 
et al. j3| discussed above and our work, tidal effects experienced by each black hole when 
close to one another. The metric is obtained by asymptotically matching a post-Newtonian 
metric for a binary system to a perturbed Schwarzschild metric modelling each black hole. 
At present, the method must be first projected onto the constraint hypersurface before nu- 
merical implementation. In addition, the method currently employs post-Newtonian results 
at the 1PN order, and considers the specific configuration of corotating BBHs. 

Turning our attention to the fourth approach, that is, construction of helical-Killing 



2a 123 and Cook 39] first showed, 



vector/quasi-equilibrium initial data, Gourgoulhon et al. 
by assuming an approximate helical-Killing vector for orbits prior to merger, a conformally 
flat spacetime and maximal slicing, that the four constraints plus one evolution equation 
simplify considerably in the CTS decomposition. The resulting set of equations is solvable 
numerically using spectral methods. Subtle problems, however, remain with the choice of 
inversion-symmetric inner boundary conditions at excision spheres around the singularities 
of the black holes - an implicit consequence of specifying a Misner topology. In this respect, 
Cook et al. |3j| recently proposed other physically-motivated quasi-equilibrium boundary 
conditions, which allow for arbitrary specifications of the conformal metric, mean curvature, 
and the shape of the excision regions. As is mentioned later, helical-Killing vector/quasi- 
equilibrium initial data sets have been studied in detail through construction of quasi-circular 



orbits and in comparison with post-Newtonian work; see [45J, |46j] and |38|,|47| respectively. 
An advantage over the double superposed Kerr-Schild method is that this approach does not 
require the specification of the extrinsic curvature. In addition, BBH dynamical evolutions 
have been performed recently for quasi-circular orbits. For example, Tl| uses a scheme based 
on excision, new gauges, the BSSN formulation [3, Q. O, 0|, and a corotating fra me, 
together with initial data constructed by means of the puncture method j3], whereas 



uses the generalized harmonic evolution scheme 12, 13, 201 based on Cook-Pfeiffer initial 

I I 

data 38 . 



In contrast to the above methods, Blanchet |1| presents an alternative set of PN moti 
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vated free data for the time-symmetric constraints, corresponding to two Schwarzschild black 
holes momentarily at rest 5 . This work proposes a particular solution to the time-symmetric 
constraints in the form of the conformal metric, 7y. The solution is such that the post- 
Newtonian expansion of the corresponding physical metric, 7^, is isometric to the standard 
post-Newtonian metric in harmonic coordinates at 2PN order as c — > +00. Furthermore, 
the solution is defined globally in space and differs from the Brill-Lindquist solution j^| for 
conformally-flat spacetimes at 2PN. It is important to note that for the numerical imple- 
mentations of BBH coalescence, the time-symmetric constraints involved relate to physically 
unrealistic events. This is because, from a general astrophysical perspective, black hole bi- 
naries prior to merger possess a velocity and, hence, impart a non-zero value to the extrinsic 
curvature to the hypersurface, i.e. 7^ 0, being the physical extrinsic curvature. 

In this respect, our work aims to establish physically realistic free data for the plunge and 
merger phase of two black holes in numerical relativity, in conformity with post-Newtonian 
results |56l ]. It follows from, and develops further, the time-symmetric approach proposed 

in El- 

Before entering into the details of our work, we remark on some of the difficulties in 
assessing how significant is the deviation arising from the assumption of a conformally flat 
metric for the merger of BBHs. Firstly, as shown by Pfeiffer et al. |42j, the resulting 
initial data sets are particularly sensitive to the choice in extrinsic curvature, as opposed 
to in the conformal 3- metric, 7^. Their work is based on a comparison of initial data 
sets based on three different conformal decompositions of the initial constraints for BBHs. 
These initial data sets have been constructed by using both the double superposition of 
two black holes 

0, y, Q 

and the Bowen-York solution j2jj. Note that although the 
superposed Kerr-Schild data sets go beyond the standard assumption of conformal flatness, 

5 The exact solutions of Misner j^, Lindquist 0| and Brill and Lindquist [^( to the time-symmetric 
constraints for conformally flat spacetimes provides an insight into the geometrostatic nature of two black 
holes, and in particular, into their topology: the solution for multiple black holes was directly inspired 
from the case of a single Schwarzschild black hole containing a Einstein-Rosen bridge (alternatively known 
as a 'worm-hole') joining two asymptotically flat universes. The Misner solution refers to two black holes, 
each containing an Einstein-Rosen bridge, which connects our universe to a second asymptotically flat 
region. In contrast, the Brill-Lindquist joins our universe to two distinct separate universes through two 
Einstein- Rosen bridges. Furthermore, it allows the computation at infinity of both the total ADM mass 
of the binary and each separate black hole, which then determines the gravitational binding energy of the 
system in the center of mass frame. 
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they disagree with post-Newtonian results at 2PN order (see Appendix |X] of this paper). 



Secondly, as is mentioned earlier, the recent results, in particular, the 



28, 



29 



38, 



38, 



ocation of the ICO 



3 



39j for co-rotating 
The agreement 



based on helical-Killing vector /quasi-equilibrium initial data 
and irrotational BBHs, show strong agreement with PN predictions 
with PN results is not surprising, considering that helical symmetry is exact with respect 
to 2PN order. Despite encouraging agreement with post-Newtonian predictions, all the 
currently used quasi-equilibrium initial data sets assume 'physically improbable' conformally 
flat spacetimes. 

C. An Overview of Our Contribution 

In this work, we propose a particular solution in the form of conformal spatial metric, 7^, 
and conformal extrinsic curvature, Kij, to the time-asymmetric constraint equations that 
represent a system of two moving black holes. The solution is chosen such that the post- 
Newtonian re-expansions of the corresponding physical metric, 7^, and physical extrinsic 
curvature, K^, are both isometric to the post-Newtonian metric in harmonic coordinates 
and post-Newtonian extrinsic curvature up to the 2PN order. Note that the post-Newtonian 
metric is formally valid in the source's "near zone" 6 . Importantly for our work, it has also 
been proved that the post-Newtonian metric arises itself from the re-expansion of a "global" 
post-Minkowskian (expansion in G) multipolar solution as c — > 00, defined everywhere in 
spacetime, including the wave zone |s7 ] . The solution for the conformal metric, 7^, is 
given, firstly, in a symmetric trace-free form, similar to the approach instigated in [l] and, 
secondly, for the first time, in Dirac coordinates. The motivation behind the choice of Dirac 
coordinates is to provide free data for constrained schema of the Einstein equations, such as 
the one based on a covariant generalized Dirac gauge 7 and spherical coordinates, as proposed 



6 The size of the npa.r zone is much smaller than the characteristic gravitational wavelength. 

7 As discussed in |58|, the advantages of implementing the covariant generalized Dirac gauge, defined as 
H2.14fl in [5^ . are numerous. Firstly, it fully specifies the coordinates in the slice X t (up to some inner 
boundary conditions for a slice containing holes); the latter property allows for the search for stationary 
solutions for the proposed set of equations, for instance, quasi-stationary initial conditions. In addition, 
the choice in gauge also results asymptotically in transverse-traceless (TT) coordinates, which are at- 
tractive for treating the problem of gravitational radiation and are analogous to the Coulomb gauge in 
electromagnetism. In addition, we show here that by relating our physical metric in Dirac coordinates, 
"fij lrac , to the post-Newtonian metric, the resulting conformal factor \]> DIRAC assumes a simple form at 
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recently by Bonazzola et al. 

Furthermore, our PN free data is constructed in such a way that it is defined at every 
point in space and differs from the Bowen-York solution at 2PN order. By demonstrating 
equivalence relationships between our PN free data and the full 2PN metric, we also prove 
that the Bowen-York extrinsic curvature is physically equivalent to the 2PN derived extrinsic 
curvature. In detail, we show that each extrinsic curvature is defined with respect to a 
different hypersurface corresponding to a distinct 3+1 foliation of the full spacetime. 

It is important to note that this work proposes only the freely specifiable parameters for 
initial data, hence, providing a free data scheme. The full physical initial data set is only 
possible after the constraint equations have been solved numerically, with appropriate inner 
and outer boundary conditions. Present numerical methods for generating a 'constraint- 
satisfying' data set from some physically motivated freely specifiable parameters include 
solving the_cpnstraints either using second-order finite difference technique s togeth er with 

J. Finite 



multigrid or successive over-relaxation j^, or using spectral methods 
regions very close to the black hole singularities are frequently removed from the computa- 
tional domain using standard excision techniques. 



D. Structure of the Paper 

The structure of the paper is as follows. Section |H] introduces the constraint equations in 
the 3+1 formalism required for presenting our approach and results. We outline the two stan- 
dard conformal decompositions of the constraints: the Extrinsic Curvature (EC) and Con- 
formal Thin Sandwich (CTS) methods, corresponding to a Hamiltonian or Lagrangian view- 
point of the problem respectively. In Section ITTT1 we first present the proposed PN free data 
in both EC and CTS decompositions, namely (7^, Al^ T , K, a) and (7^, u l \ K, N or d t K), re- 
spectively. Section ITVl subsequently outlines the proof of equivalence between our PN initial 
data and the full 2PN metric in harmonic coordinates, thus justifying the results given in 
Section imi We then briefly investigate the global properties of our proposed 'near- zone' PN 
solution and show that it is defined at every point in space. Based on our proofs, we finally 
show that the Bowen-York extrinsic curvature is physically equivalent to the 2PN derived 

2PN order; see eqn. BT7jl in Section HVBl 



S 



extrinsic curvature. Section concludes the paper with a summary of results. 



II. 3+1 CONFORMAL DECOMPOSITION 

The purpose of this section is to review the basic concepts related to the Extrinsic Cur- 
vature and Conformal Thin Sandwich decompositions of the constraints; this is necessary 
for the pres entation of the results on PN free data in Section 11111 Our account here follows 
that in |4fll . , where more details may be found. 



A. Preliminaries 



In the Cauchy formalism of Einstein's equations, a globally hyperbolic space-time is 
foliated into a set of space-like 3 dimensional hypersurfaces, X t , where each slice is a constant 
hypersurface parameterized by the time coordinate t (i 6 R). According to the standard 
3+1 decomposition, n M , a future-directed time-like unit normal to S t , is defined as = 
— iVV^t, where iV denotes the lapse function 8 . Together with the causal stability of globally 
hyperbolic spacetimes, the global time vector field, £ M , is given by t M = Nn 11 + which is 
both transverse to S t , and is normalized such that t a V a t = 1. The shift vector, thus 
trivially fulfills the condition (3 a n a = 0. Following the definition of n M , the spacetime metric 
g^ induces a purely spatial metric, 7^, defined as 7^ = g^+n^n,^, where subscripts in Latin 
and Greek characters denote three-dimensional and four-dimensional indices respectively. 

The spacetime metric components, g^, are thus expressed with respect to a preferred 
coordinate system (t,x l ) as, 

g^dx"dx u = ds 2 = -N 2 (dt) 2 + j ij (dx i + /3 i dt)(dx j + (3 j dt). (2.1) 

Together with 7^, a well-posed initial value problem must also specify how each spacelike 
hypersurface is embedded into the full spacetime. This is achieved by introducing the 

8 The lapse antidensity, a, a scalar of weight -1, is the undetermined multiplier of the scalar Hamiltonian 
constraint, which is in agreement with the canonical form of the ADM action. The relation between a 
and N is discussed later. 
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extrinsic curvature tensor, Kij, i.e., 

= -lF P n v (2.2) 

or equivalently, 

K H = -\ C nlij, (2.3) 

where C n represents the Lie derivative along the normal n 11 direction. By using the Gauss- 
Codazzi relations, the Einstein field equations are expressed in terms of six evolution equa- 
tions, 

d t K i0 = N[R %J - 2K U K] + KK l3 - torGSy + 4ttG%-(£ - p)\ - V,V,iV 

+/? Z V^- + KuVj/3 1 + KjtViP 1 , (2.4) 

and four initial constraint equations, 

R + K 2 — KijK ij = lQnGp, (2.5) 
Vj(K ij -j ij K) = 8ttGj\ (2.6) 

where eqns. (12. 5j) and (|2.6|) are referred to as the Hamiltonian and momentum constraints 
respectively. In lj'2.4jl - (|2.6jl . p denotes the matter energy density; Sy the matter stress 
tensor with S = S}, and f the matter momentum density. In addition, K = K\ gives the 
trace of the extrinsic curvature, also known as the mean curvature, and JL,- is the 3-Ricci 

n 

tensor associated with 7^. In order to avoid confusion, we follow the convention in |40J 
that covariant derivative and the Ricci tensor associated with the physical 3-metric 7^ are 
written with overbars as Vj and Rij. 

Definition ()2.3)1 additionally results in a kinematic relation between and 7^, 

dtjij = —2NKij + Vifa + Vjpi, (2.7) 

complementing the set of equations (|2.4j) - (|2.fij) . Theoretically, the constraint equations are 
thus preserved exactly under evolution (due to the presence of the Bianchi identities). 

Solutions to the most general form of the initial-value constraints for an extrinsic curva- 
ture 7^ 0, as shown by equations (|2.5p . ()2.6|) . are non-trivial, partly because the equations 
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do not fully specify 9 the dynamical, constrained or gauge nature of the components of the 
3- metric, 7^, and extrinsic curvature, Kij. Moreover, it becomes even harder to find solu- 
tions to the constraints for the specific case of two black holes. Standard procedures to solve 
the initial- value problem adopt conformal decompositions of the metric and of the extrinsic 



curvature; see ]59 



3, 



61 



6^ 



The physical 3-metric, 7^, is assumed to be conformally equivalent to a non-physical 
background metric, 7^, by a conformal factor i.e., 

Hj = (2.8) 

where \1/ is the strictly positive background (conformal) factor. 

The Hamiltonian constraint equation, f)2.5j) . as first proposed by Lichnerowicz in j^], is 
thus re-expressed for a vacuum as, 

A$ - - -m^K 2 + -^K.jK^ = 0, (2.9) 
8 8 8 

where A = V*Vj is the scalar Laplacian operator, and Vj and Rij are the covariant deriva- 
tive and Ricci tensor associated with the conformal spatial metric, 7^ (a tilde distinguishing, 
here and elsewhere, all quantities that have a conformal relationship with quantities in the 
physical space). A full conformal decomposition, however, requires introducing a confor- 
mal extrinsic curvature, K^, as developed by York [5^. Our work concerns the two most 
widely used decompositions, which are seen to be consistent with each other: a) the fully 
conformally covariant EC decomposition, introduced recently in 2^| and which improves on 
earlier non-conformally covariant decompositions, i.e. the Conformal Transverse Traceless 
(CTT) and Physical Transverse Traceless (PTT), and b) the CTS methods. Both methods 
involve expressing K%j in terms of its trace and trace-free constituents, 

= Aij + --y i:i K. (2.10) 

Each decomposition differs in its subsequent treatment of the symmetric trace-free extrinsic 
curvature, A^. The EC method is based on the 3-metric, 7^, and extrinsic curvature, K^, 
where the latter is related uniquely to the canonical momentum, 7f u . It provides a 'Hamilto- 
nian' viewpoint to the conformal decomposition of the constraints; see j4l|. In contrast, the 



Constraints l|2.5() and l|2.6|l restrict four out of the twelve degrees of freedom in (7^ , ) . 
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CTS decomposition considers the conformal metric, 7^-, and its time derivative, together 
with the mean curvature, K, and the conformal lapse function, N (where, as discussed 
in Section 111 C[ N is related to the time-derivative of K through a fifth coupled elliptic 
equation). It offers, therefore, a complementary 'Lagrangian' approach to the conformal 
decomposition of the initial data. The EC and CTS decompositions require a different set of 
free data, in particular, with respect to the extrinsic curvature or its tensor equivalent in the 
CTS instance. Both decompositions also result in, as illustrated in [4^, different physical 
initial data. Sections III Bl and III CI summarize the main results of the two decompositions. 



B. Extrinsic Curvature Decomposition - Hamiltonian formalism 

The extrinsic curvature method (see for its proposed formulation and [41] for a 
detailed review) begins by applying a weighted transverse traceless decomposition to the 



tracefree extrinsic curvature, Aij, i.e., 



A ij = -(LX) ij + Aii (2.11) 
a 

where Al£ T is a symmetric, tranverse-tracefree tensor and a is a strictly positive and bounded 
function on the 3-dimensional hypersurface, < e < a < 00 for e = constant. Given 
the vector field, X % , the notation L refers to the conformal longitudinal operator, (LX)' 3 = 
V { X j + V j X i - lY j ViX l , satisfying (LX) ij = ^- i (tX) ij . Divergence of the symmetric 
trace-free tensor, A l \ results in the following expression 10 , VjA 1 ^ = Vj[^(LX) lJ ], which can 
be solved for X % . 

Having obtained Al^ T following the substitution of X 1 in (|2.11|) . the subsequent conformal 
scaling of A l ^ T = ty~ 10 A^ T and a = \l/ 6 cr allows for, 

A ij = ^~ 10 Mgr + 4(LX)^ = ^- 10 A ij , (2.12) 



10 It involves the well-defined elliptic operator in divergence form, Vj[£r -1 (L.) y '], as discussed in |2jj and 
41 1 . For a = 1, the latter operator reduces to the vector Laplacian, (A^Y) 1 = Vj(LX) lJ . As discussed in 
23 L this is solvable on compact and on asymptotically flat manifolds, given certain asymptotic conditions 
(|59| demonstrates the existence and uniqueness of a solution to l|2.11|l for closed manifolds). In the case 
of non-compact manifolds without boundaries, boundary conditions must always be specified and their 
choice will directly affect the solution of X 1 . As remarked in [2^], there is no uniqueness property without 
boundary conditions. 
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where the weighted transverse trace-free decomposition in conformal space is A y = A^ T + 

i ( (Lx)^ 



Introduction of the weight function a in 23J ensures that the extrinsic curvature decom- 
position commutes with conformal transformations of the free data 11 . Such a specific choice 
in the conformal scaling of a also allows o to be related directly to the lapse function, N. 
As discussed in Section III CI the conformal N is a component of the free data in the CTS 
decomposition. 

In conjunction with ^j(J6~ 10 S^) = \f r ~ 10 Vj-5' u , the complete set of elliptic constraint 
equations for a vacuum are, 

Vj- (^(LX) ij ^j - ^ 6 V^ = 0, (2.13) 

A\l> - -\1>_R - — V 5 K 2 + -^AnAv = 0. (2.14) 
8 12 8 3 K J 

In this case, the set of free data comprises (7^, A^ T , K, a), which enables the solution to 
(|2.13|) and (|2.14|) for \1/ and X 1 with appropriate inner and outer boundary conditions. 
Hence, given relationships (Q, (l2~TUD . A^ T = ^- 10 i^ T , and A» = + i (jtX) ij ^j, it 
is possible to construct the physical initial data, 7^ and Kij. 



C. Conformal Thin Sandwich (CTS) Decomposition - Lagrangian formalism 
Instead of directly treating the extrinsic curvature itself (as in Section fll Bit . the CTS 

n , 

decomposition 24J considers the evolution of the metric between two neighboring hypersur- 
faces 12 . This is achieved by introducing the time-derivative of the conformal 3-metric, Uij, 

Uij = dtjij, (2.15) 

jij defined in fl2.8j) . such that = 7 1//3 <9 t (7~ 1//3 7 i: /) and j^Uij = 13 . Using a nontrivial 
conformal rescaling of both the lapse 14 , N = ^/ 6 N, and the tracefree extrinsic curvature, 

11 This is in contrast to its earlier variants (i.e. the CTT and PTT decompositions), where the conformal 
transformation and transverse-traceless decomposition are non-commutative operations. 

12 Advantageously, the method enables an understanding into the gauge choice and its subsequent evolution 
through the kinematic variables, N and (3 l . 

13 This relationship allows for the conformal metrics on both hypersurfaces to have the same determinant 
to first order in 8t 

14 This scaling is a direct consequence of the conformal invariance of the lapse antidensity a, such that 
a = a, where a and j3 l are undetermined multipliers of the constraints. When the scalar constraint is 
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A 



t3 — 2 Aij, the kinematic relation ()2.7|) simplifies to 



A ij = v[>- 10 ^' 



2N 



(2.16) 



Equation (J2.16j) implies form invariance under conformal transformations. In summary, the 
constraint equations assume the following form in the CTS decomposition, 



1 

2N 



2N 



2 a~- 

3 



i 

12 



A\I> - - ^ 5 K 2 + -^- 7 AijA ij 








(2.17) 
(2.18) 



The set of free data in this case are (7^, u %3 \ K, N). They solve for the constrained variables, 
\l/ and (3 l , in (|2.17|) and f)2.18j) . under appropriate inner and outer boundary conditions. 
Thus, the physical initial values, 7^ and K^, follow in a straightforward manner. 

We note that in practical computations, the set of four constraint equations in the CTS 
decomposition is often complemented by a fifth coupled elliptical equation relating the con- 
formal lapse, N, to the time derivative of the mean curvature, d t K. By eliminating R from 
the trace of (|2.4jl with (|2.5jl for a vacuum and after re-expressing the result in terms of their 
conformal counterparts, we obtain, 



A(A>^ 7 ) - (A># 7 

or alternatively, 

AJV + 14V* In WViN + N 
7 



8 12 8 3 



-m\d t K-(3 k d k K), 



(2.19) 



-R+-^ 4 K 2 
4 6 



4 13 



42V, In W l In * 



-V(d t K - (3 K d k K). 



(2.20) 



The free data, therefore, consists solely of pairs of variables and their corresponding ve- 
locities, (ffij,^ ,K,dtK), which is more in line with the Lagrangian approach of the CTS 
decomposition than using (fy^,^ , K, N) 41]. An additional motivation for extending the 
system of constraint equations 15 is due to the natural choice o f d+ K = in practical com- 



putations of quasi-equilibrium binary black initial data 



28, 



29 



38, 



In Section TlIIBl we 



satisfied, the ADM action results in the relationship a = Nj 1 ^ 2 , and therefore, N = ty 6 N is obtained. 
15 Despite its frequent application, uniqueness and existence proofs do not exist at present for the extended 
CTS set of equations. Pfeiffer et al. 0| recently investigated the uniqueness properties of the extended 
CTS system though without reaching any firm conclusion. Two distinct solutions for the same free data 
based on linearized quadruple gravitational waves were found in the extended system. For a given 
physical (conformally scaled) amplitude of the perturbation, the solution for the physical initial data, 7^ 
and K^, appears to be unique. 
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present the free data for both the simple and extended CTS formulations, using N and d t K 
respectively. 



D. Correspondence between the weight function a and lapse function N 

The EC and CTS formalism are equivalent to each other for the specific choice of a = 2N 
and a = 2N. Such an equivalence relationship becomes possible following the introduction 



of the weight function, a, and 



the EC decomposition; see 



3, 



w 



specifying the particular conformal scaling, o = in 



4l| for details. 



Let us assume a stationary so 



ution to Einsteins equations with timelike Killing vector 
/ such that d t gij = 0. As [23|, |4l| details, then vanishes for stationary spacetimes in 
the case a = 2N and o = 2N in the EC decomposition, where A^ T is usually identified 
with the radiative degrees of freedom. This is generally not the case for stationary, non- 
static spacetimes in the previous CTT and PTT decompositions, which do not include the 
weight function, a. Note that the standard Bowen-York free data specifies A^ T = 0; see 
27| . Alternatively, let us instead examine the particular choice a = 1, corresponding to the 
CTT decomposition 16 , in the case a = 2N and a = 2N. 

In this work, we consider both specific instances; firstly, o = 2N and a = 2N and, sec- 
ondly, a = 1 (with no assumed relationship between a and N). The second case corresponds 
to the former CTT decomposition. 



III. POST-NEWTONIAN FREE DATA 

This section presents the proposed free data, namely, (7^, A% T , K, a) and 
(jij , u l - \ K , N or d^K), based on post-Newtonian results at 2PN in EC (Hamiltonian) and 
CTS (Lagrangian) decompositions respectively. These have been introduced in Sections III Bl 
and III Cl respectively. We return in Section HVl to present our reasoning behind the choice of 
the specific form of each of the components in the above free data. 

As discussed in Section HU the EC and CTS methods offer two different perspectives in 

16 Recall that the CTT decomposition is frequently used when solving the standard Bowen-York initial data 
and is a special case within the overall EC decomposition; see Q. 
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describing how each 3-spatial hypersurface is embedded in the full space-time 17 , but do not 
differ in their definitions of the conformal spatial metric, 7^. Hence, the form of proposed 
conformal metrics, 7^, are, as expected, identical in both decompositions. 



A. EC Decomposition 

We now present each of the components in the free data (7^, A^ T , K, a). 



1. Conformal metric: jij 

Two possible options for the form of the conformal metric, 7^, are: 
(a) Symmetric- Traceless Form: 7^ 



In this case, %j assumes the form 

%G 2 m 1 m 2 d 2 g ^Gm l <4 •> AGm 2 Ki jy , „ 

= «« " ? gtfigtf + IVT". «i + "ST* ■ (3-D 



where mi and m2 refer to each of the two point particle masses respectively 18 , i.e. the 
black hole masses in our model, yi and y2 denote the black hole positions, r\ = |x— y x | 
and r 2 = |x — y 2 | represent the distances to the black holes from the field point x, and 
r i2 = |yi — ya| gives the distance between the black holes. In addition, vi = dyi/dt 
and v 2 = dy 2 /dt refer to coordinate velocities of the black holes. The term 



d 2 g 



comprises all the velocity-independent terms in (jH.ljl and represents the symmetric and 
tracefree (STF) projection 19 of the double derivative of the function g with respect to 
yi and y 2 respectively. The function g first emerged as an elementary 'kernel' for the 



post-Newtonian direct iterative works 



56 



651 ] and is defined by, 



5f(x;yi,y 2 ) = ln(ri +r 2 + ri 2 ). (3.2) 



17 In particular, the CTS method offers insights into the dynamics of the spacetime. 

18 The 'post-Newtonian' masses, mi and m2, are introduced in the post-Newtonian iteration as the coeffi- 
cients of Dirac delta functions in the Newtonian density of point-like particles. They were shown in to 
agree with the 'geometrostatic' masses associated with the Brill-Lindquist solution in the time-symmetic 
instance. 

19 i.e. Q<ij> = \{Qij + Qji) ~ \$ijQkk- 
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It satisfies the Poisson equation in a complete distributional sense such that, 

A 9 1 



(3.3) 



where A represents the standard flat-space Laplacian with respect to the field point 
x. The explicit expression for ^ is, therefore, given by, 



i9j 



d 2 9 



n\r,n 



3 



12 "12 



+ 



{n\ 2 -n\){n[ 2 + ni) 



(3.4) 



dy\dy\ r 12 (n + r 2 + r 12 ) {n + r 2 + r 12 ) 

where the notations ni = (x — yi)/Vi and n 2 = (x — y 2 )/r 2 refer to unit displacement 
vectors from x to the black holes, and ni2 = (yi — y2)A"i2 is the unit displacement 



from black hole 1 to 2. We 
Newtonian works 



56 



65 



66 



lenceforth follow the notation introduced by the post- 



63], where ^ = 



d 2 g 

dy\dy J 2 



()3.1|) may be re-expressed in a 



fully expanded form as, 



Hi 



8G 2 mim 2 



12 ""12 



rn(ri + r 2 + r l2 ] 



+ 



in 



12 



12 



n 



(n +r 2 + r 12 y 



AGm 1 



+ 



AGm 



2 „,<in,3> 



c 4 r 2 



-v o v 



(3.5) 



The precise form of (|3.1|) is chosen such that the post-Newtonian expansion (when 
c — > oo) of its corresponding physical metric 7^ is physically equ ivalent to the stan- 
dard post-Newtonian spatial metric in harmonic coordinates [56] at 2PN order mod- 
ulo a coordinate transformation. Such equivalence statements are detailed fully in 
Sections W\TK\ and ITVRl In addition, (jSHJ) is a component of the free data set 
(ffij, A l ^ T , K, a). As shown in Section llV CI this free data set refers to a solution, 
albeit approximately, to the constraints which differs from the global conformally-flat 
Bowen-York solution |2j| at second post-Newtonian order . 

(b) Dirac Coordinates 

Confining ourselves now to Dirac coordinates, the post-Newtonian derived conformal 
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metric, 7^ jrac , takes the form 20 , 



-Dirac r 



8G 2 mim 2r nTT 

— I — \m + 



AGnii 



-4 — ^1 
c 4 ri 



TT 



+ 



4Gm 2 , ,• 

—4 U 2 V 2 

c 4 r 2 



TT 



(3.7) 



where the term [ i<7j] TT denotes the transverse-traceless form of the double derivative 
of the function g (|3.2|) . with respect to yi and y2 respectively. The explicit expression 
of [ ^] TT in flUD is 



i9j\ 



TT 



(i9j) + gCijff - —Oij 



n + r 2 
ri2 



ri2 



where we refer to the notation used 
subsequently derivable, 



in 0, Q 



for D 



dy\dy. 



96 13 V 



(3.8) 

-, and the expressions 



Dg 



1 



1 



1 



0; 



d 2 g 
dy\dy{ 



2r 1 r 2 2r 1 r 12 2r 2 r 12 
1 



D 



1 1 

+ 



2rir 2 2rir i2 2r 2 ri 2 



(3.9) 
(3.10) 



Similarly, the transverse-traceless velocity-dependent terms 
2), in ()3.7j) are given by 

-I TT 



4Gmi „,j.j 



-I TT 



(and 1 



AGm 1 



— 
c 4 ri 



Gmi / 1 (^) + ^( ^(^^ + 1 Vl 2) + 3^)^) 



+n 1 J ( ^(niuij 



4 



4 



4 Vl 



(3.11) 



and 1 <-> 2, where 1 •*->• 2 denotes the exchange of particle labels 1 and 2. In more 
detail, we have chosen in ()3.7|) to fix the spatial coordinates, x l , of the conformal 3- 
metric 7^ irac on each hypersurface, St, in the generalized Dirac gauge, as introduced 



20 This may be trivially seen by applying the transverse-traceless projection tensor, TT Sfj , given by [6_ 

-\tij5 M + + VdyA- 1 - i^A-^^A- 1 ) (3.6) 



on the symmetric-tracefree conformal 3- metric, 7^ 13.111 . 
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in 0- 



The covariant generalized Dirac gauge is defined there as, 



- \ 1/3 

— I rfi Dirac 
f ' 



0. 



(3.12) 



where / and T> denote, respectively the determinant and the unique covariant deriva- 
tive with respect to a flat metric, (in an arbitrary coordinate system). Equivalently, 



(J3.12)) may be expressed in terms of the conformal metric as 
the flat metric f %3 = 5 U assumes Minkowski coordinates. 

Finally, from ()3.8|) and 1)3.11)1 . the explicit form of (J3.7)) is given as, 



d ~ij Di 



i] Dirac 



0, where 



-Dirac 



7: 



Si 



G 2 m\m 2 



()'■' 



5ri 



15 



5r 2 



+ 



12 



n.ri2 8rf 2 r 2 (n + r 2 + r 12 ) 2 



(1 



+ — + — 



n 

''"2 



+ 



12 



r 2 ri2 2rir 2 / (ri + r 2 + r 12 ) 



+n\n{ 



, (i 3) 



+ 



11 



+ 



7 



7 2 . 
— + — 

n r 12 

7 



12 



O r 2 
z ' 12 



r x r 12 %rxr\ 2 {r x + r 2 + r 12 ) 2 r 1 (r 1 + r 2 + r i2 ) 



+ ^i 2 raj 2 



(ri + r 2 + r i2 ) 2 



(ri + r 2 + r 12 ) 2 



ri 2 (ri + r 2 + r i2 ) / (ri + r 2 + r i2 



+ 



+n\n' 1 



SkJ) 



+ S lJ 



5(nit>i) 



+ 



Vl 



+ 3(riiVi)n\ v\ 



3(niWi) 2 5vx 5 



+ 1^2. 



(3.13) 



2. Extrinsic Curvature: and K 



The conformal symmetric transverse-tracefree tensor component of the extrinsic curva- 
ture, A l ^ T , is given by, 

= 0. ( 3 - 14 ) 
for both possible values of K: maximal hypersurface K = and the post-Newtonian mean 
curvature K 2 ™ derived at 2PN using results in and given by, 



K 



2PN 



Gm^nivi) Gm 2 {n 2 v 2 ) 



c 3 rf 



+ 



c 3 r 2 



+ 4 ■ 



(3.15) 



The choice of Al£ T = and = are in agreement at 2PN order with their counterparts in 
the standard Bowen-York solution, given in j^. Maximal slicing is considered advantageous 
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primarily due to its 'singularity avoiding' feature during evolution; the slicing causes the 
lapse, N, to collapse to zero in the region where a physical singularity exists; see j3| for a 
general discussion. 



3. Weight function: a 



As is mentioned in Section lHD[ we consider two specific instances; a) a = 2N and a = 2N 
and, b) a — 1. 



(a) a = 2N and a = 2N 

This concerns the class of solutions where the EC and CTS decompositions are seen to 
be equivalent. Depending on the choice in mean curvature, either K = or K = K 2 ™, 
and whether the conformal 3-metric, 7^, is specified in symmetric trace- free form or 
in Dirac coordinates (following the relationship iV = \I/ 6 iV), the conformal weight 
function based on 2PN results, <t 2PN , can take four different forms. These are: 



i) For a maximal hypersurface K = 0, <j\k=o is given as, 



a STT \ K=0 -2 = 2N 



STTi 



\K=0 



8Gm 1 1 
+ -r 



c 2 ri 



+G 2 m 1 m 2 

in a symmetric trace-free coordinate system, 
ii) Similarly, in Dirac coordinates, 

8Gm 1 



Gm i\ai \2 7 21, 35G 2 m 2 
— [3(^1*) -7v 1 ) + ^ r 

fl<->2, (3.16) 



10 



r x r 2 r 12 r 2 



~ Dirac] o n AT 

a \ K =o - 2 = 2JS 



Dirac | 



\K=0 



-2 



c 2 ri 



1 



6G mi 2 35G 2 mi 



+G 2 mim2 



( 35 



+ 



2r x r 2 r l2 r 2 



2c 4 r 2 
+ 1 <-> 2. 
(3.17) 



iii) For the 2PN derived mean curvature i^ 2PN , o"|^2pn is given as, 



<y tt |^-2pn — 2 = 2~N' 



STTi 



|^-2PN 



Gm-i 1 



c 2 ri 



+G 2 m 1 m 2 
+1^2, 



Gmi 



18 21 n 
+ + 



r\r 2 ' 2r\ 2 r 2 ' 2r 3 



2l 35G 2 m 2 



12 



2r\ 2 r 2 J 



in the symmetric trace-free coordinate system. 
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iv) Similarly, in Dirac coordinates, 

-Dirac I o r> t\t Dirac I 8Gmi 1 

(J |ft-2PN — Z = ZiV if apN — ^ = o 1 7 

c 2 ri c 4 



■[(nivi) 2 - 7v?] + 



+G 2 m m f 39 11 r i _ r i \ 
\2r 1 r 2 2r 12 r 2 2rf 2 2rf 2 r 2 / 



+ 1 «-> 2. 

(b) cr CTT in the CTT decomposition 

In this case, the conformal weight function <7 CTT = 1. 

B. CTS Decomposition 

We now present each of the components in the free data (7^, u^, K, N or d t K). 
1. Conformal metric: jij 

Due to the identical nature of the conformal decompositions of the 3-metric 7^, the 
post-Newtonian motivated conformal metric in symmetric-tracefree form, 7^, and in Dirac 
coordinates, 7^ imc , are given by (j3.1j) and ()3.7|) respectively. 



2. Time derivative of the conformal metric, Uij , and mean curvature K 

Following definition ()2.15j) and using either 7^ (j3.1j) or 7^ ()3.7j) . we propose that the 
post-Newtonian time derivative of the conformal metric, Uij, adopts the form, 

u ij = 0, (3.18) 

at 2PN for both a maximal hypersurface, K = 0, and the 2PN derived mean curvature 
K 2 ™ , given by (|3.15jl . The particular choice in Ua = and K — are in agreement with 
quasi- stationary initial conditions : 

21 Note that helical symmetry is exact with respect to Newtonian and 2PN gravity. 



28, 



29 



391. 
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3. Conformal Lapse N , or time derivative of mean curvature, dtK 

As discussed in Section Hi CI we give both the conformal lapse, N, and the time derivative 
of mean curvature, d t K, depending on whether the simple, or the extended, CTS system of 
constraint equations is to be used. 

(a) Conformal Lapse N 

Thanks to the choice in both mean curvature, K = or K 2PN , and the preferred 
spatial coordinate system of 7^ or 7^ irac , there are four different possibilities for the 
2PN based conformal lapse, N. These are: 



i) For a maximal hypersurface K — 0, N\k=o is given as, 



N 



STTi 



\K=0 



- 1 



4Gmi 1 



c 2 ri 



Gm\ ( 3 , 2 7v 2 



35G 2 m 2 
4cVf 



+G 2 mirri2 



9 



+ 



jir 2 ri 2 r 2 
in the symmetric trace-free coordinate system. 

ii) Similarly, in Dirac coordinates, 



+ l<->2, 



7\r Dirac i 

\K=0 



4Gm 1 1 



c 2 ri 



^2 f 35 
+G m x m 2 h 

V 



3Gm 1 2 35G 2 mj 
f 1 <- 



Ar\r 2 2r 12 r 2 

iii) For the 2PN derived mean curvature K 2 ™ , N\ K 2pn is given as, 



N 



STTi 



|^-2PN 



- 1 



4Gmi 1 
H 

c 2 ri c 4 



+G 2 mim2 



Griii r , . 9 9l 35(2 mf 

L 2 „ lVl ) 2 - 4«g + — 5-^1 

ri 4c 4 rf 



9 



+ 



21 



+ 



r-i 



r\r 2 4ri 2 r 2 4r? 2 ^r\ 2 r 2 



+ 1 <-> 2, 



(3.19) 



(3.20) 



(3.21) 



in the symmetric trace-free coordinate system, 
iv) Similarly, in Dirac coordinates, 

'Gmi f{n lVl ) 2 7v 2 \ 35G 2 m 2 



N 



Dirac I 



|^-2PN 



1 



AGm 1 1 



C 2 r\ 



J rG 2 m\m 2 
+1 <-> 2. 



ri 

39 11 n 
+ -: + 



(^4 ri r 2 ' 4r 12 r 2 ' 4r^ 2 ^r\ 2 r 2 



4c 4 r 2 



(3.22) 
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(b) Time derivative of mean curvature, dtK 

Since there exist two alternative slicing possibilities for our PN free data, K = or 
K = K 2PN , their time derivatives are given respectively by: i) d t K = 0, or ii) 

Gm\ , 2 2 \ G 2 mirri2 f 1 1 r\ 



2 PN Gm i for \2 2\ . G 2 m 1 m 2 ( 

d t K = "XT 3K«i - Vi ) + — s 

c 4 rj v ' c 4 ri2 V 



2r\r\ 2 2r\ 2r\ 2 r\ 



+ 1^2 + - . (3.23) 



IV. PROOF OF EQUIVALENCE WITH POST-NEWTONIAN RESULTS 

This section presents our reasoning behind the choice of free data discussed in Section ITTT1 
and, in particular, their full agreement with post-Newtonian results at 2PN order. Note that 
the free data are also chosen such that they satisfy the constraints, albeit approximately. In 
illustrating this agreement, we state first in Section flV Al our central result exhibiting the 
"near zone" behavior of the proposed solution and outline then its derivation in Section liVBI 
by considering the form of the proposed PN-derived conformal 3- metrics, 7^ and 7^* mc . In 
addition, motivations for the specific form of our PN free data are given in Section IIV CI 
by considering the lowest order perturbation of the Bowen-York solution. This allows us 
to investigate both the global structure and near-zone solution of our proposed PN based 
solution. Finally, in Section TlV D\ we show that our results imply that the physical extrin- 
sic curvature from the standard Bowen-York solution, Kf~ Y , and post-Newtonian derived 
counterpart, K 2 ™ are physically equivalent. 



A. Statement of Equivalence 

The statement of equivalence presented below is a generalization of Theorem 1 in for 
stationary black holes to moving black holes with v 1; v 2 7^ 0. 

The conformal metrics, 7^ and 7^* rac , are chosen in such a way that a post-Newtonian 
expansion (when c — ► 00) of their corresponding physical metrics, 7^- = \l/ 4 7ij and 7^ irac = 
^ -qjDirac ^^Dirac^ are physically equivalent to the standard post-Newtonian spatial metric in 
harmonic coordinates at 2PN order, that is, they differ only by a change of coordinates. 
This can be stated as, 
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A similar statement applies to Dirac coordinates. In other words, 



1 



iT c = 9tr + %tj + d £i + ° J • ( 4 - 2 ) 
In (|4.1j) and (|4.2|) . fif,,™ represents the spatial metric in harmonic coordinates truncated 
at 2PN order [56|. The remainder O [\j accounts for neglected 2.5PN and higher-order 
terms. The change in coordinates is specified by the unique spatial gauge transformation, 
or {^ Dirac ) depending on the preferred coordinate system. 

Additionally, our PN data obeys the following relationships for either maximal slicing, 
K = 0, 

7o* = SoT N + Mi + dito + o(\), (4.3) 



700 = gfiF + 2ft£o + 0{\), (4.4) 



or, mean curvature K 2 pn, 



7 0l = 9o™ + dotl + diC'o + {-s)> ( 4 - 5 ) 
7oo = gST + + o(±), (4.6) 



vC 5 / 

where £| = & + C(l/c 5 ) and g 2 ™ and g 2 f N represent the 00th and Oith component of the 

full spacetime metric in harmonic coordinates 56>j. 

The vector, which is determined by (I4.1|) - (|4.6|) . represents unique infinitesimal gauge 

transformation (i.e. x M — > x ,fl = x^+^^x"), where and are two general coordinate 

systems) at 2PN order and is explicitly given by the following components, 

Gm\ , . Gm 2 , . ( 1 \ . 

= -^r(»i*i) + ~2^M + o[^y (4.7) 

Co = + O (J;) , (4.8) 

for either mean curvature K = or K 2 pn respectively, 

G 2 m 2 o , G 2 m 2 n , ^ / 1 



for the symmetric trace-free coordinate system and, similarly, in Dirac coordinates, 
tDirac _ , _ 7G 2 m 1 m 2 n\ 3G 2 m 1 m 2 f ri + r 2 



4c 4 2c 4 (n + r 2 + n 2 ) 8c 4 V r 12 

G 2 m l m 2 (r 2 -ri\ G 2 m 1 m 2 ( r\ + r|\ 2Gm 1 , 

^ — I^^D — li" Qfc K v i r i 



c 4 V. r i2 / 48c 4 V. r-u / c 4 

4 



|j> (Afc) + ^ ( (vA0 + 1 - 2 + O . (4.10) 
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B. Relationship between the conformal metrics, 7^ and 7^"" ac , and the post 
Newtonian metric 

In order to demonstrate results (J4.1|) and ([4.2)1 . let us consider initially a post-Newtonian 
iteration of the Hamiltonian constraint, (j2.14|) and (|2.18j) . for proposed conformal metric in 
symmetric-trace-free form, 7^, at 2PN order, 

. G 2 m\m 2 „ , . Gm 1 , . Gm 2 , „ / 1 \ , , „ s 

= ( <*ffi>) + 2^?M + 2d^M + °V^ (4U) 

where we note the absence of any contributing terms from the extrinsic curvature, K^. 
As Section IIVDI shows, this is unsurprising considering that the lowest order term in the 
extrinsic curvature, K^, appears at 0(j%). Therefore, the quadratic terms of the conformal 
symmetric transverse-trace- free tensor A^ T , and mean curvature, K, occur at order 0{\) 
in (J2.14jl and ()2.18jl . The most general solution to (|4.11j) . in the full distributional sense, for 
the conformal factor, is given by, 

v 2c 4 1 3 2r 12 J 12c 4 r! v v ; ' 



Gm.2 



{3(n 2 v 2 ) 2 - v 2 2 ) + O (i) , (4.12) 



12c 4 r 2 

ip being the solution of the homogenous equation, Aip = 0. This result, when instantiated 
with v x = v 2 = 0, is consistent with the time-symmetric result of [lj], given there by 
Eqns. (2.6) and (2.7). The explicit form of ip is then obtained by comparing (j4.12|) with the 



post-Newtonian spatial metric, gff (i.e. Eqn. (7.2) in 



in the form, 



Grrii ( Grrio vi 2 \ Gm 2 ( Gm\ v 2 2 \ ,„ / 1 \ 
^ = 1 + ^( 1 -^ + 2?) + 2^( 1 -2^ + 2?) + °b)- < 413 » 

Since further 'homogeneous' terms ~ I/77 and ~ l/r 2 can be added to without affecting 
(14.11)1 . for consistency, we specify that \I/ must satisfy 1)4.11)1 in a strict distributional sense; 
see Section III A in |l|. Therefore, we do not allow the addition of such terms, here and 
henceforth, to our solutions for Poisson-type equations. 

Interestingly, by re-expressing (J4.13)) . in terms of 'Brill-Lindquist-like' constants, 07 
and a 2 , 

^1 + ^ + ^, (4.14) 
77 r 2 
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where ot\ and a 2 are determined to the relative 2PN accuracy from ()4.13|) as, 

Grri] ( Grrio v x 2 „„ ( 1 \\ , . 
01 = -W I 1 - 2^ + 2? + ° U)J 1 " 2 ' (415) 

we recover the original Brill-Lindquist conformal factor, V'bl, for the time-symmetric instance 
Vl = v 2 = at 2PN order 22 . 

In contrast, when following a similar procedure for Dirac coordinates, we find that the 
conformal factor \]/ Dimc is equivalent to the term ip, 1)4. lHj) . i.e. 

Dirac _ ^ = j + «1 + ^2 

ri r 2 

1 + gm i A _ Gm z , Xi!^ + gm2 A _ gm i , v 2 2 

2c 2 rx V 2c 2 r 12 2c 2 J 2c 2 r 2 \ 2c 2 r 12 2c 2 

+o(js)- (4-17) 

Therefore, the deviation of our post-Newtonian motivated solution from the conformally 
flat solution in the Dirac gauge manifests itself solely as a perturbation in the conformal 
metric, 7^ irac (|3.7|) . and not to the 'Brill-Lindquist '-like conformal factor, ip (j4.13)l . at the 
2PN order. 



C. Perturbation of a general global conformally flat solution 

This section considers the lowest order perturbation of the standard conformally-flat 
Bowen-York solution. It gives an insight into the specific form of the conformal 3 metric, 
7y. Note that we only consider the symmetric trace- free form of the proposed conformal 
metric, 7^ (13. lj) . This is because, as discussed in Section [IV Bl the PN based conformal 
metric in Dirac coordinates, ^f" ac (13.7)1 . incorporates fully the characteristics of the 2PN 
metric itself and does not incur a perturbative change to the Brill-Lindquist-like conformal 
factor, ip dUTJ. 



22 



tpBh is given here for completeness as, 



R-L „b-l 

V> B ~ L = 1 + — + — 

n r 2 



Gm\ / Gm 2 \ Gm 2 I ' Gmi \ ( 1 . 
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Let us introduce formally a metric perturbation, h^, to the conformal metric 7^, where 
we interpret the proposed 2PN motivated conformal metric (j3.1|) as a perturbation to the 
Bowen-York solution, i.e., 

7ij = 5 ij + h ij , (4.18) 
where h = ha = and hij is explicitly given from ()3.1j) as, 



%G mvmo d o AGm\ ... 4Gm 2 
= ~ + —r^vfvi + —r^v^vf. 4.19 



Consequently, the conformal factor, must include a perturbation, k, to the "Brill 
Lindquist-like" conformal factor, if), such that, 

* = ^ + k, (4.20) 

where r/> is given by ()4.13|) and its structure is chosen to resemble as closely as possible the 
form of the time-symmetric, Brill-Lindquist conformal factor, ipB-L (14.16)1 . 

We then assume that k, the perturbation in the conformal factor, and Kij, the full 
conformal extrinsic curvature, admit power-like expansions of the type, 



+oc / 1 \ 



n=2 



= E (^TlJ ^[n + i/ 2 ]^ (4-22) 

Subtle considerations determine the form of the above expansions. Notably, we introduce 
formally the expansions to investigate perturbations to the Bowen-York solution 27( • They 
do not, as they might otherwise misleadingly suggest, represent post-Newtonian expansions 
in the general sense. Instead, the parameter c tracks the order of the perturbation to 
the Bowen-York solution. Consequently, the lowest order of the perturbation metric for a 
non-zero extrinsic curvature, hij\Kij^o, is to be of 0(l/c 4 ). In this light, it is therefore more 
appropriate to regard the resulting perturbation equations (at low order) as successive integer 
approximations in the perturbative metric, hij, of the conformally-flat solution. Note that 
our definition ()4.22|) determines the zeroth order term in (i.e. to be of 0(1/ c 3 ), 

which is dimensionally consistent with both Kf™ and Kf~ Y , as given later by (|4.34|) and 
(14.37)1 respectively. 
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1. Analytic closed form of the linearized solution 



By considering the formal perturbation, as defined by (|4.21|) and (|4.22|) . of constraint ()2.9|) 
together with (J2.13|) and (|2.17|l . we arrive at the linearized perturbation of the Hamiltonian 
constraint, 



(4.23) 



The above is identical in form to its counterpart in the case of a linear perturbation to 
the time-symmetric Brill-Lindquist constraints, given by Eqn. (3.7) in [lj]. In contrast, the 
linear-order terms in (|4.23J) in both the perturbation metric, hij (|4.19j) . and Brill-Lindquist- 
like conformal factor, ip ()4.13j) . include terms due to the velocity of the black holes. On 
the other hand, we may identify the time-independent parts in and ip, (|4.19j) and (j4.13|) 
respectively, and refer to the results given in Section III B in [1| for solving the time- 
symmetric components in (|4.23|) . Substitution of (j4.19|) and (j4.20j) into ()4.23|) . allows us to 
determine the exact solution for Km, 



Km = —G mim2D 

+ai \ — G 2 mim 2 



g n + r 2 
6 4r 12 




+a 2 {l ^2} 



where Ai = 



dy\dy\ 



and A n 



dy\dy. 



yi and respectively, and D 



' - denote the Laplacians with respect to the source positions 
dy T dy i as before. As pointed out earlier, the results (3.8) 



and (3.11) in |l| may be used in the above. The expression for Km (|4.25|) is valid in the full 
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distributional sense, i.e. the addition of an arbitrary number of 'homogeneous' terms ~ 1/ri 
and ~ 1/Y2 is not permitted, as mentioned earlier, when solving the Poisson-type equation 
for Km. Note that k^ also tends to zero at spatial infinity (i.e. when r = |x| — > +00). 
Consistent with Eqns. (3.8) and (3.11) in Q, only terms in the first line of (|4.24j) and (J4.25|) 
contribute at the 2PN order in the conformal factor, The additional terms, which are 
proportional to 1PN constants, «i or a 2 , appear only at the 3PN order. Indeed, ()4.25|) 
reflects characteristic 3PN features, despite imposing only the initial isometric relationship 
with a 2PN expansion, (j4.1|) . In particular, the intermediate expression, ()4.24|) . contains the 
special Poisson-like solutions Hi and K±, which appear in the 3PN spatial metric, given by 



Eqn. (Ill) in 



2:-! 



. For completeness, we include here the functions, Hi and Ki, satisfying 



the following Poisson-like equations, 

AHi = 2 igjdij 



and can be found as, 



Hi 



Ax 

4 



9 



2r 12 
1 



+ D 



AKi = 2D 



o 9 



T2 



(4.26) 
(4.27) 



D 



lnri2 

n 



-D 



In?"! 
ri2 



r 2 



2r\r\ 2 



1 



2 ' 



2rfri2 2rirf 2 



Ki 



D ( — In 

J 2 



r 2 



+ 



+ 



r 2 



(4.28) 
(4.29) 



2rfr 2 2r 2 r( 2 2rfrf 2 

Finally, although ()4.24|) and ()4.25|) indicate 3PN characteristics, it is important to note 
that they do not represent complete expressions for Km at 3PN. Such an expression is only 
possible if we consider terms of 3PN order in hij (i.e. using 3PN results such as given in 



67|). which we have not attempted here. Furthermore, in order to be strictly correct up to 
and including order 0(l/c 6 ), we should instead solve for ^(3/2) by considering the 1.5 linear 
order equation to the conformally-flat Hamiltonian and momentum constraints, containing 
the first explicit appearance of the extrinsic curvature, Kij, in the hierarchy or perturbative 
equations: 



1 



lf 5K fm + l^K [3/2]ij Kii /2] 



0. 



(4.30) 



23 These occur in the non- linear potential, X. 
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0. 



(4.31) 



Such an investigation is beyond the scope of the current work. 

Finally, the fully explicit form of obtained by expanding all the derivatives in the 
result 1)4.25)1 . is given here for completeness as: 



K(i) = -G 2 mim2 



Y2r\r 2 



r\ + r 2 



12 



—Gmi 
1 



(nivi)' 
4ri 



Vl 

I2n 



Gm 



24r 12 
n 2 v 2 



4r 2 



n 



15r 2 
8r 2 r 2 2 



+Gm 2 



+ 



24rfr i2 
2r 2 





1 


r 2 J 


8r? 2 


v 2 2 " 
12r 2 _ 


+ Qti 


3 


1 



24rfr 2 4r 2 rf 2 8r 2 rf 2 



— + — 

r 2 ri 

-G 2 m 1 m 2 

2r 2 
rirf 2 



13 



12 



rfr 12 



+ 



15r 2 
Sr 2 rf 2 



+ 



2rj 



2r 2 3 

^«3^j 3 



+ Gmi 



3rir 2 r i2 

35(ni^i) 2 



(ni 2 w 2 ) 2 4(ni 2 v 2 )(niw 2 ) 



+ 



2S 2 

(^2^2) 2 _ 4(ri!t; 2 ) 2 
rxS 1 



2S 2 
+a 2 {l «-> 2}. 



S 2 

(^12f2) 2 

2r 12 S 



4(niw 2 ) 2 3(n 1 2f2)(n 2 f 2 ) 



+ 



4v 2 2 



+ 



v 2 



S 2 
2 



+ 



v 2 



35vi 2 

24r 2 
4(^1^2) (n 2 ^ 2 ) 

£2 

2v 2 2 



(7^2)^ 

2r 2 S' ' 2r 12 S' ' 2r 2 S' 3r 1 r 2 



(4.32) 



where S = (r 1 + r 2 + r 12 ). 



D. Choice of Extrinsic Curvature and Maximal Slicing 

Despite their manifestly different forms, we show here that the physical extrinsic curva- 
ture of the conformally-flat Bowen-York solution, Kf^ , and the extrinsic curvature derived 
from 2PN results, Kf™, are physically equivalent to each other modulo the infinitesimal 
coordinate transformation specified by as given by ()4.7jl - (J4. !())) ■ In particular, we show 
the following relationship, 

Kf = + d^ + o(±y (4.33) 

where £0 refers to the 0th (covariant) component of the gauge vector £ M and is specified by 
(14.7)1 . The extrinsic curvature derived from post-Newtonian results 

Q,i^f N ,at 2PN order 

is given by, 

K T = 5^ (4m««y) - + 1 - 2 + o (±) , (4.34) 
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On the other hand, the conformal Bowen-York extrinsic curvature 24 , Kfj , is given by, 

= -^[2v l(inij) - (n 1 v 1 ) - n u n y )] + 1 <-> 2. (4.37) 

We hence denote the physical Bowen-York extrinsic curvature as Kf Y = ^> 2 Kf Y , where \l/ is 
given by (|4.12j) in symmetric tracefree form and by (J4.17|) in Dirac coordinates. Therefore, 
by considering terms up to C(l/c 5 ), Kf Y = Kf Y + C(l/c 5 ). 

Note that, at first sight, the relationship ()4.33|) appears to contradict the standard trans- 
formation law of tensors. For completeness, we state the standard transformation law of 
tensors for in the instance of an infinitesimal change of coordinates, {x^} — > {x' 11 }, 
K'^ = J^J^i^p, such that the spatial components of the extrinsic curvature transform 
at linear order in are given by, 

K ij = K^ + 2d i £ j) + 0(e). (4.38) 

If we relabel = Kf™ , the difference between ()4.33j) and ()4.38j) is immediately evident. 
We stress, however, that despite this apparent contradiction, both equations are valid. The 
distinctions between the two 3-dimensional extrinsic curvatures, Kf^~ Y (which we relabel 



Strictly speaking, the Bowen-York conformal extrinsic curvature, Klj lm ~ Y , is the solution to the con- 
straints in the case of conformally flat (i.e. jij = Sij) vacuum spacetimes with a maximal hypersurface, 
K = 0, and is explicitly given by, 



r>fuiiBY 3Gmi . 

K ij = 3 2 [2v 1{i n 1:j) - {nivi) (dij - n u nij)\ 



,2 



3Ga mi 

T 2c 3 r 4 ^(xintj) + {nxVi) (S^ - n u nij)] + 1^2, (4.35) 

where we consider only the linear momentum term (with no intrinsic angular momentum) of the original 
solution. The second term in (j4.35(l corresponds to its "inversion-symmetric term" ; see [2^ . This complete 
solution i|4.35|) was chosen historically to generalize the two-sheeted topology of the Misner-Lindquist 
time-symmetric approach. As a result, the inversion-symmetric term satisfies the isometry condition for 
a field to exist on a two-sheeted manifold. The constant a in l|4.35|) denotes the radius of the inversion 
sphere, or alternatively the throat of the black hole after applying appropriate boundary conditions. 
Importantly, our approach does not concern itself directly with the topological nature of the two black 
holes. Therefore, we only refer to part of the Bowen-York solution for the most general topologies, Kfj 
(henceforth referred to as 'Bowen-York'), where 

Kf Y = ^*[2v Hinij) ~ (m«i) (fiy - n linij )] + 1 <-> 2. (4.36) 

By general topologies, we include both the two-sheeted asymptotically flat universe Misner-Lindquist 
and three-sheeted Brill-Lindquist [5^ solutions. 
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K[j in this instance) and K^, given by (|4.33jl and (|4.38jl respectively, reveal the nature of 
the 3+1 foliation of the spacetime associated with a gauge transformation 25 . 

Let us now outline the proof of the relationship (J4.33|) between the two extrinsic curva- 
tures, K[- and K' . This can be understood by first considering the infinitesimal difference, 
5n' u , between the time- like normal vectors, n' M and n' , 

8n' v = n' v - n' u , (4.39) 

where n ,fl and n ,fl denote two distinct physical vectors in the coordinate system which 
result from the transformation equations ()4.33|) and ()4.38|) respectively. From = — iVV M t 
and using the standard tensor transformation law (|4.38|) . it is possible to obtain explicit 
covariant and contravariant expressions for ()4.39|) . 

\ Nd,e + o(e) ) 

and 

*>=( -w-*F + ? + o { ?) \ 

\ Nda? + Nf3%e + Nd,e - SIMt + Ml + O(^) J 
respectively. If we then consider the quantity, SK' , the infinitesimal difference between 



Kja, and K'^, i.e. 



K» = K» + 5 K»i ( 4 - 42 ) 

where, following the definition ([2.2)1 . we find SK'^, 

5 K» = 6 HW)) 

= - (V' ( ,(K)) + (<^>' ( ,VX) + (n'")(*ny VX)) • ( 4 - 43 ) 

Substituting ()4.40|) and (|4.41|) in (j4.43J) and using 7 M „ = g^ v + n^n u , the hypersurface 
projection of the change in extrinsic curvature, SK'-, is given as, 

d(5n' ) 

M'v = — (4)r 3«)+ (4) rS(K) + ^% + ^V' k %n' 

+ n' 5n' {l Wrg.^ + n'Hn' {l ^T' k %n' , 



25 Specifically, the 3-extrinsic curvature, KL, given by the coordinate transformation 14.33JI . corresponds 
to a distinct re- formulation of the 3-foliation of spacetime. That is the 3-dimensional hypersurface T,' t , 
timelike unit vector n' a , lapse function N', and vector shift /?,', of the new coordinate system, {x ,f *}, are 
different physical entities to the corresponding E t , n a , N and [3i of the original coordinate system, {x^}. 
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where ^T'jf refers to the 4-dimensional Christoffel symbol. Finally, using n'^ = (-N', 0) 
and n' 11 = 1/N'(1, — /3 h ), we obtain at linear order in 

SK ij = - Nd ijt Wi + 9 A ~ °olij) + ^ (djjki + dij kj - dkjij) 

H ^—(djPi + diPj - d -fij) H (djjrni + da mj - d m jij) 

- ^^- (djlmi + d ilmj - d mlij ) + M^, } (-iV 2 + p k p k )) 

ffkfj cO ok am a cO 

- t—jlp-( d j)Pk) + 2N { ^ lmk + 9klm » ~ 9mlk ^ + (444) 

By considering terms up to and including 0(\), (|4.44jl simplifies to the recognizable form 
of the relationship (|4.33jl . 

^ = fy£o +°(^)- ( 4 - 45 ) 



Having proved the physical equivalence between the 2PN derived Kf™ and the standard 
Bowen-York Kf Y , it is then possible to use either the free data (Kf? N , K 2 ™) or (K^ Y , K = 
0). These results were given in both EC and CTS decompositions in Sections IlII Al and IHI Bl 
respectively. 

We finally provide complementary results to the PN free data presented in Sections IHI Al 
and IHI Bl Specifically, by assuming the widely used boundary conditions at spatial infinity, 
^jr^oo = 1 and X l \ r -> 00 = 0, we solve for the constrained variables, X 1 and j3 l , to 2PN order 
in the EC and CTS decompositions respectively: 

1. EC decomposition 

(a) Maximal hypersurface, K = 

i) & = 2N 

X 1 = -g^ (7v\ + n\(n lVl )) + 1 ~ 2 + O (I) . (4.46) 

ii) a CTT 

X* = -g^ (7v\ + nUmt*)) + 1 - 2 + O (±) . (4.47) 

(b) Mean curvature of J^ 2PN 
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i) a = 2N 



ii) a CTT 



X i = -^^v\ + l~2 + o(-\. (US) 

c 6 r\ \c° 



2. CTS decomposition 
(a) Maximal hypersurface, K = 



0* = (7«J + nKnx^)) + 1 - 2 + O (i) . (4.50) 

(b) Mean curvature with i^ 2PN 



^ = -^i«J + 1^2 + ofi). (1.51) 



V. CONCLUSION 

This work provides astrophysically realistic free data for binary black holes in numerical 
relativity, which are in agreement with 2PN results. Following the time-symmetric approach 
of lj], we propose a particular solution to the constraint equations in the form of the stan- 
dard conformal decomposition of the spatial metric and extrinsic curvature. The solution 
presented here is shown to differ from the post-Newtonian metric in harmonic coordinates 
up to 2PN order by a coordinate transformation. The solution is also shown to differ at 
2PN from the conformally-flat Bowen-York solution of the constraints, despite the singular 
nature of the proposed conformal metrics, %j and j^ irac . We recall that the post-Newtonian 
metric is not only valid in the 'near-zone' of BBHs, but also arises from the re-expansion 
of a 'global' post-Minkowskian multipole expansion when c — > oo, which is equivalent to 
a far- zone expansion when r — > oo. In addition, the post-Newtonian masses mi and rri2 
are introduced as coefficients of Dirac delta functions in the Newtonian density of point-like 
particles. Together with the formal energy and mass calculations for the time-symmetric 
instance jl| , the interpretation of our solution as a perturbation of the Bowen-York solution 
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suggests that the post-Newtonian description of the black holes as delta-function singulari- 
ties agrees with the physical masses of the Bowen-York black holes. The latter are computed 
by surface integrals at infinity and associated with Einstein-Rosen-like bridges. 

We note, however, that our solution does not include the 2.5PN term of the metric, 
associated with Newtonian radiation reaction effects. Furthermore, although our solution 
exhibits characteristic 3PN features, we do not directly use the considerably more complex 
3PN metric itself. In addition, we have only considered systems of two non-spinning black 
holes. However, spin effects are known to contribute directly to the gravitational waveform, 



and to the overall emission of energy and angular momentum of the system 
Finally, it is important to note that further studies are required on the behavior of our 
solution in the vicinity of the black holes, providing an insight into the precise nature of the 
singularity, which is necessary for practical numerical implementation of the data. 
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APPENDIX A: COMPARISON WITH THE KERR-SCHILD INITIAL DATA 



One of the most widely-used set of initial data, which assumes a deviation from a con- 
formally flat spacetime, is based on the superposition of two Kerr black holes in Kerr-Schild 



coordinates 131, 



32, 



When comparing with PN calculations, we find, however, that 
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the 'physically realistic' 2PN spatial metric, g 2PN , disagrees with the free data constructed 
by superposing two Kerr-Schild metrics together. In particular, we show that the physical 
metric, 7^ err_ child , generated from a numerically-computed conformal factor, \]> Keri '- Shlld ) j s 
inconsistent with post-Newtonian calculations at the 2PN level. 

The Kerr-Shild spacetime metric for a single black hole of mass, m, and specific angular 
momentum, a = j/m (where j is the black hole's angular momentum), is given by, 

ds 2 | K crr-Schiid = V^dx^dx u + 2H \x a )l ll l v dx il dx v , (Al) 

where rj^ is the Minkowski flat-space metric, H(x^) represents a scalar function, and 
denotes the ingoing null vector (with respect to both the background and full metric) such 
that ri^yl^ly = g^ u l^l v = (and hence, 1% = liU). For a general Kerr-Schild black hole metric 
(expressed in Kerr's original rectangular coordinates), H{x^) and Z M are, 

mr 3 . , . 

H = - A -. ttt, A2 

r 4 +(a.x) 2 ' v ' 



and 



where r is given by, 



'*Hi ™-%; + J a - x)a/r i, (A3) 



r 2 = ~ (x 2 - a 2 ) + \j\(* 2 - a 2 ) 2 + (a.x) 2 . (A4) 

From the spacetime metric (jAlj) . we obtain the 3-spatial metric, 7^, within the ADM 
decomposition, 

7 Ke rr -Schild = ^ +2W . ; (A5) 

together with the ADM gauge variables /3j = 2Hl li and a = —j===. From ()2.4|) . the 
extrinsic curvature ^Kerr-Schild j g gi ven as, 

^Kerr-Schild = _L j y . ^ + ^ _ ^Kerr-Schild] (Ag) 

The Kerr-Schild conformal metric of two black holes, 7y err_ dlild , is generated from the 
'superposition' of two Kerr-Schild coordinate systems ()A5j) . each describing a single black 
hole (see Eqn. (28) in H and Eqn. (21) in Q), i.e. 

-Kerr-Schild = + 3 ^ ^ + 3 ^ ^ (A?) 
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where the indices 1 and 2 label the two black holes. Similarly, the mean curvature is given 

as, 

^-Kerr— Schild ^^"-i^Kerr— Schild _|_ /j^-j^Kerr— Schild (A8) 

and the conformal symmetric trace-free extrinsic curvature is, 



(i ij ) Kerr " Schild = (7 fc(l ) Kcrr - Schild 



(IS fc\Kerr— Schild k Its i\Kcrr— Schild 



i / j/- fc\Kerr— Schild r k i jy- i\Kerr— Schild 



(A9) 



Where ^f^) Kerr -Schild = J_ ^fcj) Kerr-Schild , ^. + y . ^. _ ^ ^.^Kerr-Schild] &nd 

i(7 i "') Kcrr_Schlld = 2 iH(ri) ili ilj] see js2| for details. Note that for simplicity, we have 
not included the 'attenuation functions' \B and 2-B, introduced in 

For illustrative purposes, we consider here the simplest instance of two non-spinning black 
holes (i.e. a = 0), where the Kerr-Schild 3-conformal metric is given explicitly as, 

-Kerr-Schild = & + H^l^ + (A10 ) 

and the Hamiltonian constraint f)2.9j) is of the form, 

1 _ /^Kerr-SchildXS 
y\ ^jj Kerr — Schild _ ^jj>^Kerr— Schild^ _|_ \ / ^j^Kerr— Schild ^2 



(' 1 T/Kerr-Schild\-7 

■ ) (a \Kerr— Schild/ Aij\Kerr— Schild 



12 

/^Kerr— Schild ^^jj^ Kerr— Schild (All) 

where K Ke "- Schm and A* eTT ~ Schild are of order C(l/c 3 ) using {SB) and ifflfljl . 

Let us now assume that the Kerr-Schild conformal factor, ^/ Kerr - Schlld ; a t the 1PN order, 
is given by, 



^Kerr-Schild = 1 + _ + 1^ ( M2 ) 



-+(- 

c 2 \c 

which results in an explicit expression for the function Y from (jAlljl . 



Gm\ Gm 2 AGrrii BGm 2 , . 

r — — — 1 | . l Ai <jJ 

where ~ AGmi and ~ gGfm2 are additional 'homogeneous' terms, which occur when solving 
flH2j). Note that when ^Kerr-Schild = ^b-l the constants 4 _ b = 1, where V B ~ L is the 
Brill-Lindquist conformal factor given by (|4.16jl . 



Kerr-Schild 



By substituting (|A12jl and (|A13|) into the conformal relationship, 7, 

^Kcrr-Schild)4~Kerr-SctrM wg find ^ constantg A = B = 1/2 III jHHj) for the following 
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isometric relationship to be true at IPN order, 



Kerr-Schild ^1PN , r ^ i X ^ ( \ 1 A\ 

where cjd PN is the IPN spatial metric in harmonic coordinates given by Eqn. (7.2c) in 
The infinitesimal spatial gauge vector, q, is uniquely determined as, 

?i = ~^n\ - -^n\ + 0[^y (A15) 

Using A = B = 1/2 in (jA13|) and inserting v[/ Kerr -Schiid j n ^ (j2,8j) . it is immediately apparent 
that the physical metric, 7^ err_ScMLd , is not isometric with gf™ - the 2PN spatial metric in 
harmonic coordinates given by Eqn. (7.2c) in |56|. More specifically, ^Kerr-Schiid doeg 
not contain any 'interaction' terms ~ m\m2 of the two black holes. Notably, there are no 
terms involving the 'interaction' function g ()3.2|) . where we recall that g and its associated 
derivatives (such as (|3.9|) and (|3.1U|) ) are characteristic of post-Newtonian results at 2PN 
order and higher. The absence of such terms in ^-Kerr-Schild j g a d i rec t consequence from 
its construction as a linear superposition of two Kerr-Schild coordinate systems. Note that 
it is impossible to incorporate all the 2PN features of gfj into higher order terms of the 
conformal factor ^Kerr-Schild 

We, therefore, conclude that the Kerr-Schild conformal metric is incompatible with the 
inspiral physics described accurately by 2PN results. 
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